The impact of neutrino masses on the determination of dark energy properties 
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Recently, the Heidelberg-Moscow double beta decay experiment has claimed a detection for a 
neutrino mass with high significance. Here we consider the impact of this measurement on the 
determination of the dark energy equation of state. By combining the Heidelberg-Moscow result 
with the WMAP 3-years data and other cosmological datasets we constrain the equation of state 
to —1.67 < w < —1.05 at 95% c.L, ruling out a cosmological constant at more than 95% c.L. 
Interestingly enough, coupled neutrino-dark energy models may be consistent with such equation 
of state. While future data are certainly needed for a confirmation of the controversial Heildelberg- 
Moscow claim, our result shows that future laboratory searches for neutrino masses may play a 
crucial role in the determination of the dark energy properties. 



The recent results of precision cosmology and 
the measurements of Cosmic Microwave Background 
(CMB) anisotropies are in excellent agreement with 
the standard model of structure formation (see e.g. 
pj). The price-tag of this success story is a very 
puzzling consequence: the evolution of the universe 
is dominated by a mysterious form of energy, coined 
dark energy "DE" (an unclustered negative pressure 
component of the mass-energy density) with a present- 
day energy density fraction Qde — 2/3 and equation 
of state w ~ — 1 Q, 0]. This discovery may turn 
out to be one of the most important contribution to 
physics in our generation. 

Up to now, most data analysis have found a dark en- 
ergy scenario consistent with a true cosmological con- 
stant A, with w = — 1. However, the nature of dark 
matter and dark energy remains an enigma, and we 
are entering an era when their origins might be better 
understood not only through precision indirect (obser- 
vational) evidence from cosmological measurements 
but also, crucially, through direct ground based mea- 
surements of particles within and beyond the Stan- 
dard model. In this letter we discuss such a combi- 
nation of data looking at what can be learned about 
the possible interactions between neutrinos and dark 
energy from combining cosmological observations with 
the neutrino mass measurements from the Heidelberg- 
Moscow experiment (0,0, HM hereafter). 

Let us recall that in the past years mass differences 
between neutrino mass eigenstates (7711,7712,7713) have 
been measured in oscillation experiments |fj. Obser- 
vations of atmospheric neutrinos suggest a squared 
mass difference of Am 2 ~ 3 x 10~ 3 eV 2 and solar 
neutrino observations, together with results from the 
KamLAND reactor neutrino experiment, point to- 
wards Am 2 ~ 5 x 10 _5 eV 2 . While only weak con- 
straints on the absolute mass scale (Em„ = mi + 
7712 + 777,3) have been obtained from single /3-decay ex- 



periments, double beta decay searches from the HM 
experiment have reported a signal for a neutrino mass 
at > 4cr level 0, recently promoted to > 6a level by 
a pulse-shape analysis 0. As we will see in the next 
section, this claim translates in a total neutrino mass 
of Sm„ > 1.2eV at 95% c.L. While this claim is still 
considered as controversial (see e.g. 0), it should be 
noted that it comes from the most sensitive ( 76 Ge) 
detector to date and no independent experiment can, 
at the moment, falsify it. 

As very well known in the literature (see e.g. 0) 
massive neutrinos can be extremely relevant for cos- 
mology and leave key signatures in several cosmolog- 
ical data sets. More specifically, massive neutrinos 
suppress the growth of fluctuations on scales below 
the horizon when they become non relativistic. Cur- 
rent cosmological data, in the framework of a cosmo- 
logical constant, are able to indirectly constrain the 
absolute neutrino mass to Hm v < 0.75 eV at 95% c.L 
[l| and are in tension with the HM claim. However, as 
first noticed by [ljj, there is some form of anticorre- 
lation between the equation of state parameter w and 
T,m u . The cosmological bound on neutrino masses 
can therefore be relaxed by using a DE component 
with a more negative value of w than a cosmological 
constant. As we show here, the HM claim is compati- 
ble with the cosmological data only if the equation of 
state (parameterized as constant) is w < — 1 at 95%. 

While the HM claim must certainly verified by fu- 
ture experiments there are several theoretical motiva- 
tions that one should consider before discarding the 
result as simply due to unaccounted systematics. In 
the past years, indeed, the interesting idea of a possi- 
ble link between neutrinos and DE has been proposed 
The main motivation for this connection relies 
on the fact that the energy scale of DE (0(1CT 3 ) eV) 
is of the order of the neutrino mass scale. Interest- 
ingly enough, coupled neutrino-DE models predict 
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in general an effective equation of state (averaged over 
redshift) tf e // < — 1- Our result may therefore be the 
consequence of a deeper connection between DE and 
neutrino physics. 

The method we adopt is based on the publicly avail- 
able Markov Chain Monte Carlo package cosmomc 
8] with a convergence diagnostics done through the 
Gelman and Rubin statistic. We sample the follow- 
ing eight-dimensional set of cosmological parameters, 
adopting flat priors on them: the physical baryon, 
Cold Dark Matter and massive neutrinos densities, 
uJb = ^bh 2 , ijJ c = il c h 2 and £l v h 2 , the ratio of the 
sound horizon to the angular diameter distance at de- 
coupling, 8 S , the scalar spectral index n s , the overall 
normalization of the spectrum A at k = 0.05 Mpc -1 , 
the optical depth to reionization, r, and, finally, the 
DE equation of state parameter w. Furthermore, we 
consider purely adiabatic initial conditions and we im- 
pose flatness. 

We include the three-year WMAP data Q (tem- 
perature and polarization) with the routine for com- 
puting the likelihood supplied by the WMAP team. 
Together with the WMAP data we also consider the 
small-scale CMB measurements of CBI lu|, VSA 0, 
ACBAR 13] and BOOMERANG-2k2 In addi- 

tion to the CMB data, we include the constraints on 
the real-space power spectrum of galaxies from the 
SLOAN galaxy redshift survey (SDSS) JH and 2dF 
119, and the Supernovae Legacy Survey data from 
|17| . Finally, we include the Heidelberg-Moscow as in 
the recent analysis of [l8[ . 

Let us just remind that the 0^2/3 decay half-life T®" 2 
is linked to the effective Majorana mass mpp by the 
relation = m 2 /C mm T^J 2 , in the assumption that 
the 0^2/3 process proceeds only through light Majo- 
rana neutrinos and where the nuclear matrix element 
C mm needs to be theoretically evaluated. Using the 
theoretical input for C mm ( 76 Ge) from Ref. 0, the 
0^2/3 claim of @ is transformed in the 2a range 

log 10 (m/WeV) = -0-23 ± 0.14 , (1) 

i.e., 0.43 < mpp < 0.81 (at 2a, in eV). 

Considering all current oscillation data (see |Tsjl 
and under the assumption of a 3 flavor neutrino mix- 
ing the above constraint yields: 

0.0137 < n v h 2 < 0.026 (2) 

at 95% c.l. where we used the well known relation: 
n u h 2 = T,m u /93.2eV. 

Our main results are plotted in Fig.l where we show 
the constraints on the w — Hrriy plane in two cases, 
with and without the HM prior on neutrino masses. 
As we can see, without the HM prior we are able to re- 
produce the results already presented in the literature 
(see e.g. 0]), namely current cosmological data con- 
strain neutrino masses to be Hm v < 0.75 eV. However 
an interesting anti-correlation is present between the 
DE parameter w and the neutrino masses and larger 
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FIG. 1: Constraints on the w — £ plane in two cases 
with and without the Heidelberg-Moscow prior on neu- 
trino masses. 



neutrino masses are in better agreement with the data 
for more negative values of w. It is therefore clear that 
when we add the HM prior (Em„ ~ 1.8 ± 0.6 eV at 
95% c.l., again see Fig.l) the contours are shifted to- 
wards higher values of neutrino masses and towards 
lower values of id. A combined analysis of cosmologi- 
cal data with the HM priors gives —1.67 < w < —1.05 
and 0.66 < £m„ < 1.11 (in eV) at 95% c.l. excluding 
the case of the cosmological constant at more than 2a 
with Em,, = 0.85 eU, w — —1.31 and fi m = 0.35 as 
best fit values. Without the HM prior the data gives 
-1.28 < w < -0.92 and £to„ < 0.73 eV again at 95% 
c.l. with w = -1.02, Sm„ = 0.05 eV and ft m = 0.29 
as best fit values. 

The inclusion of the HM prior affects also other 
parameters. We found, at 95% c.l.: 0.916 < n s < 
0.979 (0.926 < n s < 0.989 withouth HM), 0.0209 < 
n b h 2 < 0.0235 (0.0211 < fl b h 2 < 0.0238 without 
HM), 0.302 < fl m < 0.444 (0.262 < tt m < 0.360 
without HM). It is interesting to notice that the in- 
clusion of massive neutrinos seems to further rule out 
the scale-invariant n s = 1 model. 

As we have seen the combination of a neutrino mass 
Etoj, ~ 0.85 eV and the cosmological observations in- 
dicate that the equation of state of DE w is less than 
than —1, excluding a cosmological constant. Scalar 
fields with positive kinetic energy have w > — 1 while 
phantom fields can have w < — 1 but they have 
a negative kinetic energy and many fundamental the- 
oretical problems. Here, we prefer to follow the ap- 
proach of interacting DE [7j • Interacting DE are mod- 
els where the dark energy interacts with other parti- 
cles, as for example dark matter or neutrinos. The 
net effect of this interaction is to change the apparent 
equation of state of DE 0. An observer that sup- 
poses that the DE has no interaction sees a different 
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evolution of DE as an observer that takes into ac- 
count for the interaction of DE. This effect allows to 
have an apparent equation of state w < — 1 for the 
"non-interaction" DE observer even though the true 
equation of state of the DE is larger than -1. 

To sees this more clearly let us define the energy 
density and pressure of the scalar field as p^ — ^<j) 2 + 

V(<t>)> P<j> — \y — V{4>) and the equation of state 
parameter = p^/ p<t>, where the potential V(</>) does 
not include the interaction with dark matter. If there 
is no interaction between DE and dark matter, w<p 
gives the complete evolution of DE and w > —1 . We 
now include an interaction term between dark matter 
(or neutrinos) with <j> via the function /(</>) which gives 
an interacting dark matter energy density |2f)j | 
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where f = /(</>„) and a Q — 1 at present time. In this 
case dark matter no longer redshifts as a~ 3 since the 
evolution of f((f>) will also contribute to the redshift. 
The evolution of piu and <j> are given by 

Pim + 3H(piM + pdm) = BlMlL. 4> (4) 

O- Jo 

4> + 3Hj> + V = -^-t (5) 
a d So 

where the prime denotes derivative w.r.t. 0, i.e. 
V = dV/dcj), f = df/d(f>. The total dark matter does 
not need to coincide with pjm ■ This would be the case 
if we want to interpret pim as the energy density of 
neutrinos since we know that they cannot give the to- 
tal amount of dark matter. However, since neutrinos 
are massive they certainly contribute to dark matter. 

It was pointed out in Q that the apparent equation 
of state, i.e. the equation of state of DE if we had 
assume that there was no interaction, is 



W0 



u ap 



1 - x P<t>a A \ } 



In this case the noninteracting DE and dark matter 
observer sees a standard evolution, p m = ~'SHp m and 
Pde = —3H pde^X + w a P )- We see from eq.© that 
for / < fo we have x > and w ap < , which allows 
to have a w ap less than -1. 

The effect of an apparent w ap less than -1 has a 
stronger effect on small redshifts when the DE dom- 
inates. This effect is measured by the SNIa and the 
actual values for the redshifts of these supernovae are 
mostly in the range < z < 1.2. So, let us expand 
the function f(cf>(a)) as a function of the scale factor 
around a = 1, 



SW =So + 



dS 



da 



LO-i) + 



(7) 



For generality and presentation purposes we assume 
that the scalar field is already tracking, i.e. we take 



u>0 constant, and then the energy density is given by 

P4> = P4>oa- 3{1+w ^ = (2/(1 - u>*))V($, where we 
have used that the kinetic energy can be expressed as 
Ek = (1 + w$)/(l — Wj,)V. Taking the derivative of 
p^ w.r.t. a we can relate d(f>/da to the potential V, 
d(f)/da = — 3(l+w t / > )V/(aV'). Using this expression of 
d<fi/da, p^ as a function of a and eq.JJI), we can write 
x in eq.©, in terms of redshift z=l/a-l, as 



= 3A 



1 



- W4, 
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+ 



3w<b 



(8) 



where we have defined A = —CIimo^Im/^ with 
Ajm = Sol So, A0 = V„/V . A positive x requires 
A > 0. Notice that the evolution of w ap and x depends 
on z only via the term z/(l + z) 3w <t> in eq.© and once 
and fl^o are fixed the value of w ap is determined 
by the present day values of VL IMo , So,V ,Sq, V only 
through A. 

Since w ap is a function of z it is better to use 
the weighted average equation of state to compare 
the models with the observational data. The average 
equation of state is defined by 



/ W a ptt DE dz 

J Q DE dz 



(9) 



where the integral runs from z = to z — 1.2. The 
effective w e ff is then a function of Clfo, w<f, and A. 

We show in Fig.© the evolution of func- 
tion of z (dashed line) for fl,p = 0.65, A — 0.35 (i.e. 
fl IMo = 0.35 if A /M /A = -1), = -0.98. We see 
that w ap decreases with increasing z and becomes less 
than -1 at z = 0.3 and it is w ap — —1.6 at z = 1.2. 
We also show in Fig.@ the behavior of w e ff as a 
function of A (solid line) with the same parameters 
fi^o = 0.65, Wff, — —0.98. With increasing values of 
A, w e ff becomes more negative and for A — 0.1 we 
find w e ft = — 1 and at A = 0.84 we have w e ff = — 1.3 
as required by the cosmological plus HM data. Finally, 
if we assume that the interacting matter is only due 
to neutrinos with the total amount of neutrinos today 
given by the central values of the CMB plus HM anal- 
ysis Y>m v = 0.85 eV then n IMo h 2 = Sl v h 2 = 0.009 
and Ajm/A^ = —40 for w e ff = —1.3. 

We have seen in a model independent study that 
using interacting DE it is possible to obtain w e f / less 
than -1, consistent with the values given by the cosmo- 
logical data plus HM. Future high -z baryon acoustic 
oscillation surveys, in tandem with high-z supernovae 
surveys should provide a powerful mechanism to look 
for such deviations from w = — 1. 

To summarize and conclude, we have considered in 
this letter the cosmological implications of the con- 
troversial Heidelberg-Moscow result. We have found 
that a scenario based on a cosmological constant is 
unable to provide a good fit to current data when a 
massive neutrino component as large as suggested by 
HM is included in the analysis. A better fit to the 
data is obtained when the DE component is described 
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FIG. 2: We show function of z (dashed line) and 

w e ff as a function of A = —Q.im \im/^4> (solid line) with 
integration limits < z < 1.2 for f^ = 0.65, W4, = —0.98. 

with an equation of state w ~ —1.3, with w; < — 1 at 
more than 95% c.l.. As far as we know, this is the 
only dataset able to exclude a cosmological constant 
at such high significance. 

There exists, therefore, a significant tension be- 
tween the indirect, observational measurements lead- 
ing to the LCDM scenario and the direct HM obser- 
vations. Rather than implying one should rule out 
evidence from the direct measurements purely on the 
basis of disparity with the indirect observations, this 
tension suggests we should keep our minds open to 
alternative dark energy scenarios beyond a cosmolog- 
ical constant. This, together with the fact that the 
energy scale of DE (0(1O~ 3 ) eV) is of the order of 
the neutrino mass scale, may suggest for a link be- 
tween neutrino physics and DE that must certainly 
be further investigated. Systematics can be present 
in the HM data and a more conservative treatment 
(see 0) would lead to a better agreement with a cos- 
mological constant. However, phantom models with 
w < — 1 would still provide a better fit to the data. 



On the other hand, using a more conservative ap- 
proach towards cosmology, by, for example, combin- 
ing HM only with the CMB dataset, would provide 
even larger values of £?7i„ and more negative values 
for w. Recent combined analysis with Lyman-a for- 
est data (0,01) hnply tight constraints on neutrino 
masses (£m„ < 0.2eV), seemingly at discord with 
the HM result, and also in some tension with CMB 
data alone. Future larger scale Lyman-a surveys and 
refinements in the analysis, addressing systematic un- 
certainties and sensitivity to modeling assumptions, 
will allow a better assessment of how these tensions 
will be resolved. 

If the HM result is correct, a signal in the range 
m v ~ few x 10 _1 eV is clearly expected, and could be 
found in the next-generation Karksruhe Tritium Neu- 
trino Experiment (KATRIN) [U, which should take 
data in the next decade, with an estimated sensitivity 
down to ~ 0.2 eV. Cosmological data will also reach 
a 2a accuracy of about ~ 0.1 eV on Em„ in the fu- 
ture [fj mostly thanks to Planck satellite CMB experi- 
ment expected to launch in 2008, and future precision 
high-z supcrnovae and baryon oscillation dark energy 
surveys. 

A determination of the absolute neutrino mass scale 
will therefore not only bring relevant information for 
neutrino physics but may be extremely important in 
the determination of the dark energy properties and in 
shedding light on a possible neutrino-dark energy con- 
nection. Future direct particle detection and indirect 
astronomical experiments will scrutinize this interest- 
ing hypothesis. 
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